Numerical renormalizat ion- group study of the Bose-Fermi Kondo model 
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We extend the numerical renormalization-group method to Bose-Fermi Kondo models (BFKMs), 
describing a local moment coupled to a conduction band and a dissipative bosonic bath. We apply 
the method to the Ising-symmetry BFKM with a bosonic bath spectral function rj(u>) oc ui s , of 
interest in connection with heavy-fermion criticality. For < s < 1, an interacting critical point, 
characterized by hyperscaling of exponents and w/T-scaling, describes a quantum phase transition 
between Kondo-screened and localized phases. Connection is made to other results for the BFKM 
and the spin-boson model. 
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I 



I 

C 

O 

o 



(N 
> 

o 

o 
in 
o 



o 
o 



X 



Dynamical competition between local and spatially ex- 
tended degrees of freedom provides a possible mechanism 
for the intriguing non-Fermi liquid behaviors observed 
in cuprate high-temperature superconductors [l| and in 
heavy-fermion systems near a quantum phase transition 
(QPT) 2]. The essence of this competition is embod- 
ied in a new class of quantum impurity problems: Bose- 
Fermi impurity models of a local moment coupled both 
to a fermionic band and a dissipative bosonic bath. 

The key physics underlying heavy-fermion QPTs — 
interplay between screening of /-shell moments by 
conduction electrons and magnetic ordering of those 
moments due to the Ruderman-Kittel-Kasuya-Yosida 
interaction — is contained in the Kondo lattice model 
(KLM) H H 0. The KLM has been studied via an 
extended dynamical mean- field theory (EDMFT) H @ 
that maps the lattice onto a self-consistcntly determined 
BFKM in which bosons represent the fluctuating effective 
magnetic field generated by other / moments. Assuming 
the bosonic bath has a sub-Ohmic spectrum t](uj) oc lu s 
(s < 1), a perturbative (in e = 1 — s) renormalization- 
group study has found two types of antiferromagnetic 
ordering transition: a QPT of the usual spin-density- 
wave type 0; and a locally critical QPT, characterized 
by s = + , at which divergence of the spatial correla- 
tion length coincides with critical local-moment fluctua- 
tions. There is growing evidence for the latter scenario 
in several materials, particularly from neutron scatter- 
ing results on CeCus.gAuo.i |8|. Similar mappings to 
an effective Bose-Fermi impurity model have been ap- 
plied to the disordered KLM 9] and to a t-J model of 
a highly incoherent metal close to a Mott transition, of 
possible relevance to the pseudogap phase of the cuprates 
jlfij . The BFKM also describes certain dissipative meso- 
scopic qubit devices, e.g., a noisy quantum dot, where the 
bosonic bath represents gate- voltage fluctuations [llj . 

A satisfactory description of all these problems hinges 
on proper solution of the BFKM. The identification of 
a locally critical QPT in [§| was based on perturbative 
results for e — ► 1. Most other techniques developed for 
conventional (pure-fermionic) impurity models [l2j are 



inapplicable or suffer from significant limitations. For 
example, nonperturbative quantum Monte Carlo stud- 
ies [a] of the KLM have reached conflicting conclusions 
concerning the existence of a locally critical QPT due to 
inherent limitations in accessing the lowest temperatures. 

To meet the need for solutions of the BFKM that are 
reliable for all bath spectra and all temperatures T, we 
have turned to Wilson's numerical renormalization-group 
(NRG) technique 0], which has hitherto provided a con- 
trolled treatment of a range of fermionic impurity and 
lattice problems 0] > an d was recently applied (lRflil to 
the spin-boson model (SBM) of a two-level system cou- 
pled to a dissipative bosonic bath Below we report 
a nontrivial development of the NRG that incorporates 
coupling of an impurity to both a conduction band and a 
bosonic bath. We give results for the BFKM with Ising- 
symmetry bosonic couplings, the case most relevant to 
CeCus.gAuo.i and a noisy quantum dot. For bosonic 
bath exponents < s < 1, the model displays criti- 
cal properties identical to those of the SBM, with expo- 
nents satisfying hyperscaling relations and w/T-scaling 
in the dynamics. We show directly the destruction of the 
Kondo resonance at the quantum critical point (QCP). 
Importantly, this work opens the way for future decisive 
studies of the KLM, obtained by imposing EDMFT self- 
consistency on the BFKM solution. 

We study the Hamiltonian H — Hp + Hb + Hi nt , with 

Hf = ^2e k cl a c ka , H B = ^2 W q4>\^q, 

k.cr q 

H mt = JS-s c + gS z J2^i + ^- q ), (1) 
q 

which describes a spin-^ local moment S interacting 
with both the on-site spin s c of a fermionic band, and, 
via S z , with a bosonic bath. The essential physics of 
the model, based on perturbative RG and large-iV re- 
sults 0, 0, 113, H3, is the competition between the 
Kondo coupling J, which tends to screen the impurity 
moment, and the bosonic coupling g, which favors an 
unscreened spin. The fermionic band is assumed to be 
featureless with a density of states p(e) = po for |e| < D. 
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The impurity-boson interaction is embodied entirely in 
B(ui) = Trg 2 r](u}) — Bq uj s for < oj < luq, where 
s > — 1. In the EDMFT problem of future interest, s 
will be determined self-consistently from the form of the 
local spin susceptibility xi c(w) and will in general take 
a non-Ohmic value (s^ 1) y|. To connect with previous 
work |20|. we write -Bo = {K^g) 2 , and for convenience we 
set u>q = D = 1 in the following. 

NRG scheme. The band (bath) continuum of energies 
\e\ < 1 (0 < cj < 1) is replaced by a discrete set, ±A~™ 
(A - ™) for n = 0, 1, 2, where A > 1 parameterizes the 
discretization. Eq. may then be cast in a chain form, 
in close analogy to the pure-fermionic NRG, introducing 
a discretization error that vanishes as A — > 1 0] : 



H = 



y ] ^nfnufna + T n(fn a fn-l,a + fn—l,afnv) 
n,cr 

J2 [ E mblb m + T m (blb m -! + fcLl&m)] (2) 



P° J y^/oo-^'/oa' • s + K o9\ — S z {b\ + b 
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Here, foa (bo) annihilates the unique combination of 
fermions (bosons) that couples to the impurity. The 
tight-binding coefficients £„, r„, E m and T m (m,n — 
0, 1,2, ...) encoding all information about the band and 
the bath are determined via a Lanczos procedure. The 
decay of e n and r„ as A~™/ 2 for large n is what allows the 
iterative solution of the conventional Kondo problem via 
diagonalization of progressively longer chains . Since 
the bath spectral function r](uj) vanishes for uj < the 
coefficients E m and T m decay faster, as A _m for large to, 
as discussed previously for the SBM 15]. 

We solve the discretized BFKM via an iterative scheme 
that adds a site to the fermionic chain at each step, 
but extends the bosonic chain only at every second step. 
Since T?, n ~ r n , this ensures that fermions and bosons 
of the same energy scale are treated at the same step. 
The number of particles per bosonic site is in principle 
unlimited, but is restricted in practice to a finite maxi- 
mum Nb- For the conventional Kondo model, the basis 
must be truncated after only a few iterations, saving only 
the lowest N s eigenstates of one step to form the basis (of 
4N S states) for the next step. Truncation is required even 
sooner in the BFKM, where the basis grows by a factor 
of 4:(Nb + 1) at each step at which bosons are added. 

It is not obvious a priori that the approach outlined 
above should capture the physics of the BFKM QCP , 
but it is validated by the results presented below. These 
results are converged with respect to the truncation pa- 
rameters, 8 < N b < 12 and 500 < N s < 2000 typically 
sufficing. Calculated critical exponents prove insensitive 
to the choice of A (see Table UJ, so to reduce computer 
time A = 9 was used, except where noted otherwise. 




FIG. 1: (a) Schematic phase diagram for the BFKM, showing 
Kondo (K), localized (L), and quantum-critical (QC) regimes 
separated by a crossover scale T* oc \g — g c \" that vanishes 
at the QCP. (b) and (c) Correlation length exponent v as a 
function of bath exponent s. Crosses (at s — 0.2 and 0.8) show 
representative results for the SBM. Dashed lines are linear fits 
discussed in the text. 



Critical coupling. For s < 1 the NRG flow of effective 
couplings has two stable limits: a Kondo fixed point at 
J = oo, g = 0, describing a phase in which conduction 
electrons screen the local moment and the single-particle 
spectrum exhibits a Kondo resonance; and a "localized" 
fixed point at J = 0, g = oo, where the impurity dynam- 
ics are controlled by the coupling to the bath. A third, 
unstable fixed point — the QCP at J* ^ 0, g* ^ — is 
reached for bare couplings lying precisely on the bound- 
ary g = g c (J) between the Kondo and localized phases. 
For s > 1, by contrast, the Kondo fixed point is reached 
for all J > 0. At s = 1 the critical point is Kosterlitz- 
Thouless-like ^jji while for — 1 < s < it is nonin- 
teracting. We restrict attention henceforth to the range 
< s < 1 over which the QCP exhibits properties (shown 
below) characteristic of an interacting critical point. (For 
the case s = of interest in connection with local criti- 
cality, logarithms replace power laws at the QCP.) 

Correlation length exponent. The crossover from the 
quantum-critical regime to a stable low-temperature 
regime defines a scale that vanishes at the critical cou- 
pling: T* ~ \g — g c \ u [Fig. ^a)]. The correlation length 
exponent ^(s), which can be determined from the NRG 
many-body spectrum or from properties such as the local 
susceptibility or the single-particle spectrum, diverges for 
,s — > + and s — ► 1~, indicative of qualitative changes at 
these limiting values of the bath exponent. For small s, 
v ~ 1/s [fit in Fig. njb)], which becomes asymptotically 
exact as s — ■> 0, while for s — * 1 the exponent approaches 
v = 1/^/2(1 - s) + C [fit in Fig. dc)]. Our v values are 
in essentially exact agreement with NRG for the SBM 
[HI and, for s -> 0, with perturbative RG for the BFKM 
expanding about s = |23| . 

Local magnetic response. Central to the EDMFT of 
future interest is the local susceptibility xioc^ the re- 
sponse to a magnetic field h acting solely on the impurity. 
Fig. shows the static local susceptibility x\oc(T;uj = 0) 
vs T for ,s = 0.2. The two phases are readily distin- 
guished. For g < g c , screening of the impurity moment 
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FIG. 2: Static local susceptibility xioc(T;uJ — 0) vs T for 
s — 0.2, and poJ — 0.5, where Kog c ~ 0.355. The inset shows 
the vanishing of Txioc(0;0) for g ^ gt ■ 



is signaled for T < T* by a constant value xioc(T;lj = 
0) oc (g c — g) 1 that diverges as g — > g~ . For g > g c , 
by contrast, there is free-moment behavior with a Curie 
constant Txioc(T;co — 0) oc (g — g c ) X (see inset to Fig. 01. 
At criticality, it is found that 

Xio C (T;u = 0,g = g c )<xT- x , (3) 

with a; = s for all < s < 1, as predicted by the e- 
expansion [5nj . We have also calculated exponents de- 
fined by M Xoc (g > g c , T = 0, h = 0) oc (g - gt)' 3 and 
Moc(5 = g c ,T = 0) oc |/i|V<5, w here M loc = (&) serves 
as an order parameter for the localized phase. Table |l] 
shows exponents obtained for two representative cases. 

Starting from a scaling ansatz for the critical part of 
the free energy, F crit = Tf(\g - gd/T 1 /" ', \h\/T b ), it is 
readily shown that 5 = (l+x)/(l — x), A = 2/3 = u(l — x), 
and 7 = vx. The exponents obtained numerically for 
< s < 1 obey these hyperscaling relations, indicating 
that the critical point is interacting |]j. 

We now turn to the imaginary part of the dynamical 
local susceptibility, x[oc( w )- Figure. El shows xlod^'T — 
0, g = g c ) vs oj for s — 0.2. The exponent y, defined via 

xL(b>; T = 0, g = g c ) cx sgn(«), (4) 

is found to satisfy y = s for all < s < 1. The static 
local susceptibility x\oc(T) is also shown, from which it is 
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0.2 4.99(5) 


0.200(1) 


0.665(5) 


1.99(2) 


1.00(1) 


0.8 2.11(2) 


0.803(3) 


0.106(2) 


0.209(1) 


1.67(2) 



TABLE I: Critical exponents for A = 9, s = 0.2 and 0.8. 
Parentheses surround the estimated nonsystematic error in 
the last digit. The insensitivity to A is illustrated by the 
s = 0.2 values v = 4.90(5), x = 0.200(1) for A = 5, and 
v = 4.9(1), x = 0.199(4) for A = 3. 



evident that x = y, which likewise holds over the entire 
s range. Due to truncation the NRG method is known 
to be unreliable for |ct>| < T [2^|. However, the result 
y = x for < s < 1 (obtained in for the particular 
case s = + ) is consistent with cj/T-scaling at the QCP 
and reflects the interacting nature of the fixed point 0). 
Specifically, wc find for lo <C Tk that 

T KX L(u,T;g = g c )~ ^ (|) , (5) 

with <j) s a universal function of its argument as demon- 
strated in the lower part of Fig. A similar result has 
recently been obtained for the multi-channel BFKM in 
the large- N limit [i^. 

An important observation is that all critical exponents 
of the Ising-symmetry BFKM are in numerical agreement 
with those of the corresponding sub-Ohmic SBM (l5lll6|. 
This clearly indicates that that the two QCPs belong to 
the same universality class [as implied in Ref. 6(a)]. 

Spectral function. Finally, we calculate A{ui) = 
— —\va,{{d;d')) UJ for the Ising-symmetry Bose-Fermi An- 
derson model described by 

H mt = X^ £d + l Un d-cr) n da +V^2(dlc ka +h.C.) 
o k.a 

+ \g{nd] -n dl )^{<t) q +0t g ). (6) 
q 

Here, = d\,d a is the spin-cr occupancy of an impu- 
rity level of energy e^, U is the on-site interaction, and 
V is the hybridization between the impurity and the 
conduction band. The coupling to the bosons is again 
through the z component of the effective impurity spin. 
In the regime Z7, je^l 3> r = np V 2 where charge fluc- 
tuations are suppressed, this Anderson model maps onto 
the BFKM under a Schrieffer- Wolff transformation • 
The QCP is manifest in the single-particle spectrum as 
a collapse of the central Kondo resonance. Fig. 21 shows 
A(u>) vs lo close to and at the QCP. On the scale of the 
Hubbard satellite bands [Fig. Ufa)], the spectra are barely 
distinguishable. Fig.^Jb) shows the same spectra on the 
scale of the central feature. For B < B c the integrity 
of the Kondo resonance is preserved; it remains pinned 
at the Fermi level, A{u — 0) = (7rr ) _1 , as required by 
Fermi-liquid theory for the conventional Anderson model 
flij l . Close to the transition, the lowest-frequency (Fermi- 
liquid) behavior 1 — 7tT q A(uj) oc lu 2 crosses over above a 
scale \w\ ~ oj* to a power law characteristic of the QCP it- 
self. For B Q = S ,c [dashed line in Fig.QJb)], uj* = 0, i.e., 
the spectrum remains pinned but the lowest-frequency 
behavior is non-Fermi-liquid in nature. For Bq > Bq c 
the Kondo resonance collapses and A(lj — 0) decreases 
with increasing Bq. Near the critical coupling there is a 
crossover, seen as a weak maximum in A(to) at \oj\ = uj* 
[Fig. life)], to the quantum-critical behavior described 
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FIG. 3: Upper-right axes: Xi'oc(u>',T = 0) (dashed) and 
Xioc(T; u = 0) (dotted) at the QCP for s = 0.2 and p J = 0.5, 
showing that x — y — s; see Eqs. @ an d Q- Lower-left 
axes: Finite-T scaling of Xioc{w,T) at the QCP as described 
by Eq. @. Data are shown for s — 0.2, poJ = 0.1 (thin line) 
and 0.25 (thick line) with g — g c in each case, and for tem- 
peratures T/T K = 10" 1 , 10" 2 , 10" 3 , 10~ 4 , 10" 5 , 10~ 10 and 
10 -20 . The curves have a common form for uj /T <C Tk/T. 




FIG. 4: Single-particle spectrum vs u> for bath exponent 
s = 0.7, impurity parameters U = — 2sd = 0.002 [such that 
A(lu) — A(—lj)], hybridization strength To = 1.5 x 10~ 4 , and 
discretization A = 3. Shown for Bo — So,c = (dashed 
line), ±10 -4 on (a) the scale of U and (b) the scale set by 
the central feature, (c) Shown on a logarithmic scale for 
B Q - B , c = (dashed line), ±10~ 5 , ±10" 4 . 



above. The correlation length exponent may also be cal- 
culated from the vanishing of uj* at the transition, and 
agrees with the results of Fig.^to within numerical error. 

Summary. We have developed the NRG method to 
treat coupling of a quantum impurity to both a conduc- 
tion band and a bosonic bath, as described by Bose-Fermi 
Kondo and Anderson models. Our approach provides 
a good account of the critical properties of the Ising- 
symmetry BFKM with a power-law bath spectrum. Over 
the range of bath exponents < s < 1, the competition 
between fermionic and bosonic couplings gives rise to an 
interacting fixed-point, exhibiting hyperscaling of criti- 
cal exponents and ui/T scaling in the impurity dynam- 
ics, which belongs to the same universality class as the 
critical point of the SBM. The method can readily be ex- 



tended to models having multiple bosonic baths, such as 
the XY and isotropic BFKMs. This work also opens the 
way for EDMFT treatments of the Kondo lattice down to 
T = that will shed much-needed light on the possibility 
of local criticality in heavy- fermion systems. 

We thank Q. Si for useful comments. This work was 
supported in part by NSF Grant DMR-0312939. 
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